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Engineering  Department 

ABSTRACT 

A  theoretical  framework  for  modeling  the  characteristics  of 
railgun  plasma  armatures  is  developed  taking  into  account  non-ideal 
transport  properties  and  equation  of  state.  The  electrical  resistivity 
of  the  plasma  is  calculated  as  a  sum  of  contributions  from  electron-ion 
and  electron-neutral  collisions.  Non-ideal  correction  to  the  ionization 
potential  and  pressure  in  Saha’s  equation  and  the  equation  of  state  is 
included.  For  electron-ion  resistivity,  we  use  the  basic  formulation  of 
Khalfaoui  for  highly  non-ideal  plasma,  smoothed  and  calibrated  with 
respect  to  the  results  of  Kurilenko  &  Valuev,  Vorobiov  et  al  and 
Kovitya.  For  electron-neutral  resistivity,  we  use  a  formulation  which 
allows  the  fitting  of  an  infinite  half-power  series  to  experimentally 
determined  energy  spectrum  of  electron-neutral  collision  frequencies  for 
a  given  atom. 

Next,  a  steady-state  magneto-gasdynam i c  model  improved  over 
earlier  works  is  developed.  In  this  model,  the  net  Lorentz  force  exerted 
on  the  plasma  armature  is  exactly  accounted  for. 

Finally,  we  apply  the  model  to  study  the  radiation  emanating 
from  the  plasma  armature.  It  is  found,  on  the  one  hand,  that  the 
radiation  from  the  plasma  armature  under  the  conditions  of  most  railgun 
experiments  conducted  in  the  past  is  sufficiently  high  to  cause 
ablation.  On  the  other  hand,  the  modelling  results  also  indicate  that 
there  indeed  exist  domains  of  operating  parameters  in  which  radiative 
ablation  of  the  wall  materials  can  be  completely  eliminated.  Whether 
these  favorable  domains  of  operating  parameters  can  be  realized  in 
practice  remains  a  development  of  the  future,  though  preliminary  experi¬ 
ments  conducted  in  a  companion  program  sponsored  by  DoE  and  experiments 


I ntroduct i on 


1  . 

Theoretical  models  for  raiigun  arcs  have  been  introduced  by 
McNab[l],  Powell  and  Batteh[2],  Th i o  [3— 5^ ,  Powell  and  Batteh[6], 

Powel I [7]  and  Batteh [8] .  These  studies  have  been  useful  in  planning 
plasma  driven  raiigun  experiments  in  recent  years.  They  have  been 
particularly  useful  in  designing  suitable  raiigun  diagnostics  and  in  the 
analysis  of  experimental  data. 

In  these  earlier  works,  liberal  use  was  made  of  simple 
expressions  for  the  plasma  transport  properties  and  the  equation  of 
state,  based  essentially  upon  ideal  plasma  behavior.  In  these  studies, 
the  plasma  armatures  were  estimated  to  have  rather  high  temperatures 
(30,000  K  -  60,000  K)  and  relatively  low  pressures  (100  MPa  -  200  MPa). 
At  these  temperatures  and  pressures,  ideal  plasma  behavior  was  a 
reasonable  assumption,  making  these  studies  internally  consistent. 

Through  the  early  experimental  work  of  Hawke  et  a  I  [9] ,  Thio  et 
a  I [ 1 0] ,  Bedford  et  al[ll],  Clark  and  Bedford[12]  and  Bedford[13],  it 
soon  became  clear  that  ablation  of  the  wall  materials  was  commonplace  in 
railguns  in  which  the  arc  was  generated  by  exploding  metallic  foils  of 
relatively  low  mass.  This  ablation  adds  considerable  mass  to  the  arc  in 
a  relatively  short  time,  whence  the  arc  begins  to  expand  and  cool 
leading  to  the  formation  of  a  relatively  cool  arc  at  high  pressure.  The 
arc  evolves  and  enters  the  regime  of  non-ideal  plasmas.  Parker[14] 
later  emphasized  the  important  role  this  ablation  played  in  degrading 
raiigun  performance. 

The  first  objective  of  the  work  reported  here  is  to  develop  the 
appropriate  theoretical  capability  so  that  we  can  extend  our  investiga¬ 
tion  of  the  behavior  of  ra i Igun  arcs  into  the  realm  of  cold,  high 


pressure  arcs.  The  theoretical  construct  we  shall  develop  in  this  paper 
is  based  upon  non-ideal  plasma  theory. 

The  second  objective  is  to  introduce  a  refinement  on  the 
momentum  equation  and  boundary  conditions  used  in  the  earlier  works. 

With  this  refinement  the  accelerating  Lorentz  force  can  now  be 
accurately  accounted  for. 

The  third  objective  is  to  apply  the  model  developed  to  study  the 
radiation  emanating  from  the  plasma  armature  over  a  broad  range  of 
operating  parameters  in  a  railgun.  Before  the  advent  of  this  work,  there 
were  many  conjectures  on  the  question  of  what  are  the  dominating 
physical  processes  responsible  for  ablation  of  the  wall  materials  in 
railguns.  Many  of  these  speculations  remain  to-day  and  need  to  be 
explored.  In  this  regard,  what  has  been  accomplished  from  this  work  is 
that  radiation  has  been  theoret i ca I  1 y  established  as  an  important 
physical  process  to  consider  in  the  context  of  wall  ablation  in 
railguns.  This  does  not  preclude  the  possibility  of  other  mechanisms 
being  present  and  possibly  being  important  or  even  more  important.  What 
it  implies  is  that,  if  wall  ablation  in  a  railgun  is  to  be  eliminated, 
radiative  ablation  is  a  significant  component  which  needs  to  be 
e I i m i nated . 


2.  A  MAGNETOGAS DYNAMIC  MODEL 


Powell  and  Batteh  [2 , 6-8]  and  Thio[3-5]  have  previously  developed 
steady  state  or  quasi-steady  state  1-D  models  for  the  plasma  armature. 
The  model  we  shall  develop  in  the  present  study  shares  a  similar 
starting  point  with  these  earlier  models  but  differs  from  them  in 
several  essential  details.  Investigations  on  the  structure  of  railgun 
arcs  have  also  been  undertaken  by  Stainsby  and  Bedford [38]  and 
Ma rsha I  I [39 , 40] ,  in  which  the  f  ine  structure  of  the  armature  is 
emphasized,  especially  in  the  boundary  layer  near  the  wall.  Our  study 
here  focuses  on  the  interior  of  the  arc  (the  arc  column)  sufficiently 
far  away  from  the  wall  so  that  the  arc  may  be  reasonably  considered  as  a 
continuum.  The  boundary  layer  between  the  wall  and  the  arc  column  is 
not  considered  here 

Figure  1  shows  the  geometry  of  the  model.  The  axes  are  chosen  so 
that  the  z-axis  is  in  the  direction  of  projectile  travel,  the  y-axis  is 
normal  to  the  rails,  and  the  x-axis  is  horizontal  into  the  paper.  We 
assume  that  all  physical  variables  can  vary  at  most  in  the  z-direction 
(1-D)  and  that  the  entire  system  comprising  the  plasma  armature  and  the 
projectile  has  reached  a  steady  state  and  moves  at  a  constant 
acceleration  driven  by  a  constant  total  current.  The  momentum  equation 
for  the  plasma  armature  in  that  case  can  be  written  as: 


dg 

oz 


c  a  =  jB 


(2) 


where  p  is  the  gas  pressure,  o  is  the  gas  density,  a  is  the  acceleration 
of  the  entire  system,  j  is  the  current  density  and  B  is  the  magnetic 
field  Ampere’s  law  takes  the  simple  form, 
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Figure  1  —  Geometry  of  the  plasma  model. 

where  a  (4ir  x  10  ^H/m)  is  the  permeability  of  free  space.  A'e  require  an 
equation  of  state  to  relate  the  pressure  to  the  density.  The  ideal  gas 
equation  of  state  was  used  in  the  earlier  models  of  T hi o ,  Powell  and 
Batteh  Herein  we  introduce  the  non-ideal  gas  pressure  as  a 
perturbation  on  the  ideal  gas  pressure,  and  write  the  equation  of  state 
in  the  form, 

x  o  kR  T 

p  =  {1  *  .1  a.Z.  *  g(7/Z)>  {<) 

P 

where  g ( 7 / Z )  1 s  a  slowly  varying  function  of  the  non-ideal  parameter 

1/2  (see  expression  (18)),  and  represents  the  non-;deal  oerturbation  to 
the  true  pressure,  2  is  the  ionization  fraction  of  the  ion  in  the  i-th 


and  cannot  be  accurately  approximated  by  1-D  analysis  without  some  , 

modifications  to  take  account  of  the  3-D  effects.  In  particular,  the 
field  at  the  forward  boundary  of  the  plasma  armature  is  not  zero. 

These  reasonings  led  Thio[3,4]  originally  to  suggest  an  alter¬ 
native  approach.  While  accepting  the  fact  that  Maxwell’s  equations  in 
1-D  require  that  the  boundary  magnetic  field  should  equal  its  value  at 
infinity,  Thio  decided  to  discard  the  requirement  that  the  field  at  the 
forward  boundary  B  must  be  zero.  This  opens  up  precisely  the  needed 
degree  of  freedom  necessary  to  obtain  a  better  estimate  for  the  unknown 
magnetic  field  at  the  rear  boundary.  In  Thio[4],  the  rear  boundary 
magnetic  field  was  obtained  by  actually  averaging  the  field  across  the 
mid-line  of  the  bore  in  between  the  rails  produced  by  the  current 
flowing  in  the  finite-height  rails  and  in  the  plasma.  The  rails  were 
assumed  to  be  infinitely  long  and  the  ra i I  current  was  assumed  to  flow 
only  on  the  inner  rail  surface.  The  method  gives  a  much  better  estimate 
of  the  accelerating  force,  typically  within  20%  of  the  true  value. 

Batteh[8]  later  improves  on  this  model  by  averaging  the  field  across  the 
bore  of  the  gun,  not  only  for  the  field  at  the  rear  boundary,  but  also 
for  each  arc  cross-section  (x-y  plane)  normal  to  the  z-direction.  In 
doing  so,  Batteh  obtained  a  much  improved  momentum  equation  for  the  arc, 
but  sacrificed  much  of  the  simplicity  of  the  1-D  approach.  The  resultant 
equation  was  an  i ntegro-d i f f erent i a  I  equation. 

In  this  paper,  we  suggest  yet  another  approach  which  is  designed 
to  give  exactly  the  correct  value  of  the  Lorentz  force  without 
sacrificing  the  mathematical  simplicity  of  the  1-D  approach,  and  at  the 
same  time  retaining  strict  consistency  with  Maxwell’s  equations  in  1-D. 

We  begin  with  the  expression  for  the  Lorentz  force  given 
globally  from  field  theoretic  consideration  in  the  form, 

F,  (t)  =  \  L ’ (t)  I?(t)  ♦  M’  (t)  I  (t)  I  (t) 


(11) 


where  Ir(t)  is  the  current  flowing  in  the  inner  rails  and  I  is  the 
current  flowing  in  the  augmentation  windings,  is  the  inductance 

gradient  of  the  inner  rails,  M’^  is  the  mutual  inductance  gradient 
between  the  inner  rails  and  the  augmentation  windings.  The  values  of 
the  inductance  gradients  L*r  and  pertaining  to  a  specific 

application  can  be  accurately  obtained  either  by  an  actual  measurement 
with  an  LCR  meter  on  a  given  gun  or  by  computation  using  a  2-D  or  3-D  EM 
field  analysis  code  (Patch  et  al[15],  Thio  et  al[16]).  The  magnetic 
fields  at  the  rear  and  forward  boundary  are  now  chosen  to  satisfy  (8), 

(9)  and  the  exact  value  of  the  Lorentz  force  (11).  The  results  are, 
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(13) 


where, 

L;ff(t)  =  l;  ♦  2  [ia(t)/ir(t)]  m>2  (14) 

is  an  effective  inductance  gradient  for  the  system  comprising  the  inner 
rails  and  the  augmentation  winding. 

The  resultant  model  not  only  is  capable  of  including  effects  of 
flux  augmentation  (series  or  trans-augmentation),  it  is  also  ideally 
suitable  for  modeling  non-steady-state  behavior  of  plasma  armatures 
subject  to  time-varying  current  and  time-varying  rail  inductance 
gradient.  For  this  purpose,  time  dependence  has  been  explicitly 
indicated  in  the  above  expressions. 

We  need  an  equation  for  the  temperature  T  of  the  plasma 
armature.  The  temperature  is  determined  by  the  radiative  equilibrium  of 


state  of  ionization,  kR  is  the  Boltzmann  constant  (1.38  x  10  J/deg  K) 

B  -27 

and  (Dp  is  the  mass  of  the  atomic  species  (A.W.  x  1.66  x  10  kg).  A 

brief  discussion  on  the  non-ideal  perturbation  term  in  this  equation  is 

given  later  in  Section  4. 

For  a  uniform  current  and  temperature  distribution  within  the 
arc,  the  above  set  of  equations  can  be  solved  analytically,  provided  we 
treat  the  non-ideal  function  g  as  a  constant  parameter.  The  solution 
yields  the  pressure  distribution  within  the  arc  to  be, 


p  (Z)  =  2  {b  (1  -  e~kla^)  -  k  af} 

(kja) 

where, 


B  (Za)  kj. 


*iZi  *  g(7/Z)>  V 


and  £  =  z  -  zg,  zg  is  the  z-coor dinate  of  the  rear  boundary  of  the 
plasma  armature.  Using  the  equation  of  state  (4),  the  density 
distribution  can  be  determined  from  (5),  and  the  total  mass  for  the 
plasma  armature  can  be  obtained  as, 


m,  =  (^)  (r2)  T'1  I  I  2 

A  vah'  '■bg  p 


where , 


b2  =  {abi  [1  +  kTIT  (e  1  '  1)]  -  kial} 

V k  ^  a  1  y  1 


b_  =  fl  +  E  aZ.  + 


In  both  the  expressions  for  the  pressure  distribution  and  the  total 
armature  mass,  the  magnetic  field  B  =  B(z  )  at  the  rear  boundary  enters 
as  a  parameter  and  needs  to  be  determined. 

Equation  (3)  expressing  Ampere’s  law  requires  that  the  fields  at 
the  rear  and  the  forward  boundaries  be  connected  by  the  relationship, 

B(zb)  -  B(za)  =  ^  (8) 

where  1^  is  the  total  current  flowing  through  the  plasma  armature,  and  h 
is  the  bore  height  (the  height  of  that  portion  of  the  rails  from  which 
the  arc  collects  its  current).  Also,  the  net  Lorentz  force  exerting  on 
the  plasma  armature  as  a  whole  is  rigorously  given  by  the  volume 
integral  of  j  *  B,  which  from  the  1-D  Maxwell’s  equations  gives 

FL  =  k  K  -  Ba>  (9) 

In  Powell  and  Batteh[2],  the  forward  boundary  magnetic  field  was  set  to 
zero.  The  rear  boundary  magnetic  field  B^  follows  from  expression  (8). 
Upon  substituting  this  boundary  magnetic  field  into  expression  (9), 
Powell  and  Batteh[2]  obtained  for  the  net  Lorentz  force, 


F  =  -  (^)  I2 
L  2  ;h  '  r 


(10) 


While  the  procedure  can  be  shown  to  be  rigorously  correct  for  the  case 
of  a  true  1-D  geometry  in  which  the  rails  are  infinitely  high,  it  tends 
to  overestimate  the  net  acceleration,  typically  by  a  factor  of  2  or 
more.  As  can  be  seen  from  (10),  the  equivalent  rail  inductance  gradient 
L’  is  u  (1.257  /iH/m)  for  the  case  of  a  square-bore  railgun  (w=h)  , 
whereas  in  practice  the  actual  value  of  L’  for  a  square-bore  railgun  is 
typically  03  -  04  jH/m 

The  basic  cause  of  this  discrepancy  is  that  the  magnetic  field 
near  the  boundaries  of  the  plasma  armature  is  essentially  3-D  in  nature 


the  plasma  armature.  The  ohmic  dissipation  in  the  plasma  armature  is 
balanced  exactly  by  the  radiative  loss  from  the  arc  to  the  side  walls 
and  through  the  front  and  the  rear  boundaries, 
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where  R^  is  the  net  resistance  of  the  plasma  armature,  and  we  have 

allowed  an  emissivity  £^  for  radiation  to  the  rails,  a  separate 

em’rssivity  £■  for  radiation  to  the  insulator,  and  an  emissivity  £  for 
i  ■'a 

radiation  through  the  rear  boundary,  and  another  emissivity  through 
the  forward  boundary.  The  plasma  resistance  is  calculated  from  the 
express i on , 


(17) 


where  I  is  the  length  of  the  plasma  armature,  h  is  the  height  of  the 
bore,  w  is  the  width  of  the  bore  (separation  between  the  rails),  and  rj 
is  the  mean  resistivity  of  the  plasma  armature. 

To  make  further  progress,  we  need  expressions  which  would  give 
reasonable  estimates  of  the  plasma  resistivity  and  emissivity. 


3.  RESISTIVITY  OF  RAILGUN  PLASMA  ARMATURES 


In  the  search  for  optimum  plasma  armatures  to  avoid  the  ablation 
problems  of  the  early  rail  guns,  we  were  led  gradually  but  steadily  to 
the  realm  of  plasmas  which  exhibit  strongly  non-ideal  behavior. 

For  most  of  the  plasmas  considered  in  the  present  study, 
significant  ionization  is  present.  The  degree  of  ionization  is 
typically  in  excess  of  205?.  In  these  plasmas,  the  deviation  from  ideal 
behavior  arises  mainly  from  the  collective  electrostatic  interactions 
among  the  charged  and  neutral  particles,  in  opposition  to  the 
randomizing  effect  of  their  thermal  energy.  The  ratio  of  the  mean 
coulombic  energy  of  interaction  of  the  particles  to  their  mean  thermal 
energy , 
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where  e  is  the  electronic  charge,  2  is  the  effective  charge  of  the  ions, 
and  n^  is  the  density  of  all  charged  particles,  is  a  good  indicator  of 
the  severity  of  the  non- idea  I  behavior.  Note  that  we  have  used  the 
density  of  all  charged  particles  including  electrons  in  the  definition 
of  7  rather  than  the  more  common  practice  of  using  only  the  electron 
density  ng.  The  higher  the  value  of  the  non-ideal  parameter  7,  the 
greater  is  the  deviation  from  ideal  behavior.  Generally,  ideal  plasma 
theory  is  applicable  only  when  7  <<  0.1.  For  values  of  7  >  0.1,  non¬ 
ideal  plasma  theory  is  required.  At  a  fixed  temperature,  the  degree  of 
non-ideal  behavior  increases  as  the  electron  density  increases.  For 
electron  density  below  the  critical  density, 
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where  m  is  the  electron  mass  (9.11  *  10  kg)  and  h  is  the  Planck 
6  -34 

constant  (6.624  x  10  ns),  the  electron  energy  distribution  in  the 
plasma  remains  essentially  classical  in  accordance  with  the  Maxwellian 
energy  distribution.  Beyond  the  critical  density,  the  electron  energy 
distribution  becomes  non-Maxwe I  I i an  and  is  governed  by  Fermi  statistics 
with  the  gas  then  becoming  a  quantum  plasma[18]. 

In  this  study,  our  computations  are  limited  to  classical  non¬ 
ideal  plasmas.  Most  of  the  plasmas  we  examined  have  density  at  least  a 
factor  of  5  below  the  critical  density,  though  some  come  as  close  as 
within  a  factor  of  2  of  the  critical  density,  when  the  validity  of  the 
classical  Maxwellian  distribution  becomes  questionable.  We  have  not 
extended  our  search  for  plasma  armatures  into  the  domain  of  quantum 
plasmas.  There  is,  however,  no  apriori  reason  for  not  doing  so. 

Currently  an  active  area  of  research,  the  work  on  the 
characterization  of  non-ideal  plasmas  has  produced  voluminous 
collections  of  largely  unsystematized  experimental  data  and  a  great 
variety  of  theoretical  approaches.  We  make  no  claim  to  the  complete 
coverage  of  available  literature  and  results  on  the  subject.  A  good 
review  of  the  subject  of  electrical  transport  properties  for  these 
plasmas  can  be  found  in  the  works  of  Kha I f aou i  [17]  and  others  [18-21] . 

The  electrical  resistivity  f]j  of  a  classical  plasma  is  the 
result  of  all  the  collisions  among  all  the  particles  present,  and  can  be 
given  in  terms  of  the  collision  frequencies  or  equivalently  the 
relaxation  times  of  the  various  species  of  particles  by, 


where  1/  and  v  .  are  the  collision  frequencies  of  electrons  with 
en  e  i 

neutral  particles,  and  ions  respecti vely,  and  their  mean  values 
indicated  by  the  brackets  <  >  has  been  taken  with  respect  to  the 
normalized  Maxwellian  energy  distribution, 

f  (E)  =2  *_1/2  (kT) "3/2  E1/2  exp  (-E/kgT)  (21) 

3.1  COULOMBIC  RESISTIVITY 

For  the  contribution  to  the  resistivity  from  the  interaction  of 
charged  particles  with  charged  particles,  we  based  our  calculation 
mainly  on  the  theory  of  Kha I f aou i  [18]  who  gives  a  comprehensive  model 
covering  both  the  non- idea  I  classical  and  quantum  plasmas.  For  much  of 
the  domain  of  our  computation,  we  find  that  the  Khalfaoui  resistivity 
shows  only  a  weak  dependence  on  the  non-ideal  parameter  7 /Z,  and  can  be 
reasonably  approximated  by  the  following  simple  expression, 

Ve]  =  3.5  x  io2  ST  exp  (-27/32)  T"3/2  (23) 

showing  clearly  the  predominant  3/2-  power  dependence  on  the 
temperature . 

Despite  the  3/2-power  law  of  the  temperature,  the  difference 
between  expression  (23)  and  Spitzer’s  resistivity  in  functional  form  and 
numerical  values  is  fundamental  and  significant.  Spitzer’s  resistivity 
involves  a  logarithmic  dependence  on  the  ratio  of  the  Debye  length  Ln  to 
the  impact  parameter  p^  where  Tp^  plays  the  role  of  a  scattering  cross- 
section.  This  arises  from  setting  the  Debye  length  as  the  cut-off 
distance  for  coulombic  interaction  in  approximating  the  relaxation  times 
of  electron-ion  encounters  in  deriving  Spitzer’s  res i sv. i v i ty .  The 
approximation  makes  sense  so  long  as  there  is  a  sufficiently  great 
number  of  particles  within  the  Debye  sphere  (radius  Lp)  to  shield  the 
nuclear  charge  of  the  ion  from  outside  the  Debye  sphere,  so  that  nearly 
all  the  electron-ion  encounters  take  place  within  the  Debye  sphere. 
However,  at  sufficiently  high  electron  density,  due  to  the  long  range 


nature  of  the  electrostatic  interaction,  a  significantly  large  amount  of 
coulombic  scattering  occurs  outside  the  Debye  sphere.  This  is  reflected 
numerically  in  the  value  of  the  scattering  radius  Pg  being  similar  in 
magnitude  or  greater  than  the  Debye  radius  (A  <  1)  .  The  Debye  length 
can  no  longer  serve  as  the  cut-off  distance  in  evaluating  the  relaxation 
time  integral.  Consequently,  Spitzer’s  expression  for  the  resistivity 
fails. 

In  the  range  of  particle  density  and  temperature  relevant  to  our 
present  study,  we  have  found  that  Spitzer’s  expression  can  underestimate 
the  coulombic  resistivity  by  as  much  as  a  factor  of  5  compared  with 
expression  (23). 

3.2  THE  ELECTON-NEUTRAL  SCATTERING  RESISTIVITY 

As  the  plasmas  under  consideration  are  not  fully  ionized,  a 

significant  number  of  collisions  could  occur  between  the  electrons  and 

the  neutral  particles.  Their  contribution  f?  to  the  overall 

'en 

resistivity  should  therefore  be  appropriately  taken  into  account, 
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The  principal  difficulty  here  is  the  evaluation  of  the  mean 
e I ectron-neutra  I  collision  frequency.  The  theoretical  problem  consists 
of  solving  the  Schrod i nger ’ s  equation  for  the  scattering  of  electrons  by 
a  many-e I ectron  atom,  calculating  the  phase  shift  and  evaluating  the 
scattering  cross-section  using  the  standard  phase  shift  formula[33].  We 
shall  not  follow  that  procedure  here.  Instead,  we  shall  make  use  of 
empirical  expressions  for  the  dependence  of  the  collisional  frequency  on 
electron  energy  deduced  from  actual  experimental  measurements. 

For  the  gaseous  atoms  of  interest,  we  have  found  it  possible  to 
approximate  the  energy  dependence  of  the  collisional  frequency  between 
the  electrons  and  atoms  within  experimental  limits  of  scatter  and 


discrepancy  using  the  following  functional  form: 
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where  the  b’.s  see  coefficients  chosen  to  best  fit  the  experimental 

data,  E  is  the  electron  energy  in  J  and  n  is  the  number  density  of  the 

-3  n 

neutral  particles  present  in  m  Taking  the  average  with  the 
Maxwellian  energy  distribution  (21),  the  contribution  to  the  plasma 
resistivity  from  the  electron-atom  scattering  is  computed  in  the  form, 

f]  =  g  r3/2  (26 
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where  the  function  q  is  given  by, 
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Using  the  data  from  Frost[23]  and  Cambel[24],  a  table  of  values  for  the 
various  coefficients  and  atoms  have  been  compiled  as  shown  in  Table  1 
and  used  in  this  study. 

Before  we  can  use  the  results  of  this  section,  the  degree  of 
ionization  which  determines  the  densities  of  the  charged  species  at  a 
given  temperature  and  pressure  needs  to  be  determined.  For  this 
purpose,  we  need  an  equation  of  state  and  an  expression  giving  the 
effective  ionization  potential  for  the  atomic  species  under  non-ideal 
cond i t i ons . 


/  / 


4.  THE  NON-IDEAL  EQUATION  OF  STATE 


The  equation  of  state  for  high  density  matter  has  been  a  subject 
of  extensive  research  for  the  last  fifty  years.  Computational  methods 
have  grown  steadily  from  relatively  simple  models  to  highly  sophisti¬ 
cated,  complicated  ones.  The  starting  point  in  most  of  these  methods  is 
a  model  for  the  potential  of  the  atom  in  a  given  ionization  state. 
Through  this  potential,  charged  and  neutral  particles  collectively  or 
individually  interact.  The  Helmholtz  free  energy  associated  with  this 
interaction  is  derived,  from  which  the  equation  of  state  is  determined. 
Recent  work  includes  those  of  L  i  berman  [25]  ,  Kerley  &.  Abda  I  I  ah  [26]  ,  and 
Rinker[27]  etc.  More  historical  work  includes  those  of  Cowan  4 
Ashkin[28]  using  the  Thomas-Ferm i -D i rac  model  of  the  atom.  Extensive 
tabulations  of  the  equation  of  state  for  a  great  number  of  elements  over 
a  wide  range  of  pressures  and  temperatures  can  be  obtained  from  the 
library  of  SESAME[29]. 

A  related  effect  is  the  lowering  of  the  ionization  energy,  a 
subject  which  by  itself  has  attracted  many  investigators  and  is  richly 
endowed  with  I i terature [30-34]  . 

The  detailed,  accurate  quantum  dynamical  computations  of  the 
equation  of  state  is  beyond  the  scope  of  this  study,  nor  is  it  required. 
In  the  present  context,  the  pertinent  effect  due  to  non-ideal  behavior 
of  the  plasma  is  a  small  reduction  in  the  gas  pressure  and  ionization 
energy  due  to  the  attractive  electrostatic  interaction  between  the 
charged  particles.  A  reasonable  estimate  for  these  corrections  could  be 
obtained  using  classical  or  sen i -c I  ass i ca I  methods. 

In  the  realm  of  classical  or  semi -c I  ass i ca I  methods,  a  variety 
of  theoretical  results  are  available.  Thus,  for  a  weakly  non-ideal 


plasma  this  pressure  reduction  has  been  calculated  within  the  framework 
of  the  Debye-Huckel  theory  (see  for  example  Landau  k  Lifshitz[34]  or 
Zel’Dovich  k  Raizer[35]  p.  216).  That  theory  requires  that  there  are 
sufficient  number  of  particles  within  the  Debye  sphere  to  justify  the 
statistical  treatment  on  which  the  Debye-Huckel  theory  is  based.  As  the 
electron  density  increases,  the  mean  number  of  particles  within  the 
Debye  sphere  decreases  to  below  unity,  at  which  point  the  Debye-Huckel 
theory  completely  loses  its  physical  meaning.  This  corresponds 
approximately  to  values  of  7  >  0  1 .  Unsold  in  1948[30]  introduced  the 
so-called  1-particle  theory,  which  was  analysed  later  by  Ecker  k 
Kroll[31]  to  make  physical  sense  only  for  very  high  density,  typically 
for  7  >>  1.  Even  in  that  range,  Ecker  k  Kroll  had  found  that  Unsold’s 
result  could  give  unrealistically  high  estimates  for  the  coulombic 
interaction  energy,  which  is  reflected  in  an  excessive  lowering  of  the 
ionization  potential.  Ecker  and  Kroll  [31]  attempted  a  smoothing  of  the 
results  of  these  two  diverse  approaches.  A  unification  of  these 
classical  theories  based  on  a  refinement  of  the  theory  of  Ecker  and 
Kroll  is  recently  given  by  Thio[36].  For  our  immediate  purpose,  we 
shall  adopt  the  model  formulated  in  Thio[36],  which  gives  simultaneously 
an  expression  for  the  pressure  correction  and  the  lowering  of  the 
ionization  energy . 

T  h i o  [36]  gives  the  pressure  correction  for  a  non-ideal  plasma 
over  the  full  range  of  the  non-ideal  parameter  7  from  0  to  approximately 
1  as, 
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where  Lp.  is  the  Debye  radius  which  is  gi*en  by, 


and  f“o  is  the  mean  interparticle  separation  for  the  charged  particles, 
defined  as: 
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In  the  above  expressions,  n.  denotes  the  number  density  of  the  i  th  ionic 
species,  and  Z.  the  ionic  charge  in  units  of  electronic  charge  (Z-  =  i). 

Combining  the  pressure  correction  with  the  ideal  gas  pressure, 
we  may  write  the  equation  of  state  as, 
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where  the  function  g (7/ Z)  is  given  by, 

9(7/2)  =  '[a!  (7/2) 3/2  [1  -  a2  (7/2)1/2]-1/3 

with, 
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The  associated  lowering  of  the  ionization  potential  for  the 
transformation  of  the  m-th  ion  to  the  (m+l)-th  ion  is  given  by  Th i o [36] 
as , 
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where  the  function  f  is  given  by 


f(S)  =  1  -  *  2?2  In  (1  ♦  |) 

Saha’s  equation  for  the  (m+l)-th  state  of  ionization  is: 
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where  U  is  the  partition  function  of  the  ion  in  the  mth  state  of 
m 

ionization. 

The  results  developed  in  this  section  can  now  be  combined  with 
those  of  Section  3  to  calculate  resistivities  of  real  plasmas  at  a  given 
temperature  and  pressure,  which  can  be  compared  with  actual  exp  rimental 
measurements.  In  Figure  2  we  show  this  comparison  for  high-pressure 
lithium  plasma  for  which  experimental  data  have  been  given  by 
Vorobiov[22]  for  its  resistivity  as  a  function  of  temperature  (up  to 
15000  K)  and  pressure  up  to  100  MPa.  The  agreement  with  the 
experimental  values  are  reasonable  considering  the  uncertainty  in  the 
experimental  data.  In  the  same  figure,  we  have  also  shown  the 
calculated  resistivity  for  a  H  plasma  as  a  function  of  temperature  at  a 
specified  pressure  of  73  MPa,  and  compared  it  with  the  computat i ona I 
result  of  Kovitya[37],  Again  the  agreement  between  the  two  computed 
results  are  reasonable  considering  the  lack  of  relevant  experimental 
data  in  this  case . 
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Figure  2  —  Resistivity  of  Li  and  H  plasmas, 
present  theory  and  others 


Comparisons  of  results  o 


5.  THE  MODEL  SOLUTIONS  AND  SCALING  RELATIONSHIPS 


The  main  results  from  the  formal  solution  to  the  magnetogas- 
dynamic  equations  can  now  be  assembled  as  follows: 


Dynamical  equilibrium: 

:  1  =  (2h/ fiw)  (b3/b2)  T  I"2  mA 

(33) 

Radiative  equilibrium: 

:  l\  =  2  as  g2  1  T4 

(34) 

Dissipative  element: 

ra  =  w  yhi 

(35) 

Electrical  transport: 
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and  I,  m^,  T,  R^  and  f]  are  the  length,  mass,  temperature,  resistance, 
and  resistivity  of  the  plasma  armature  respectively.  1^  is  the  current 
flowing  through  it  and  the  inner  rails  while  I  is  the  current  in  the 
augmentation  windings.  L*  is  the  inner  rails  self  inductance  gradient 
and  's  the  mutual  inductance  between  the  inner  rails  and  the 

augmentation  winding.  We  have  written  the  relationships  (33)  -  (36)  in 
such  a  way  as  to  expose  the  main  functional  relationship  among  the 
principal  variables:  the  plasma  armature  length  I,  mass  m^,  temperature 
T,  resistance  R^,  resistivity  rj ^  and  the  plasma  current  1^,  where  the 
functional  dependence  on  these  variables  is  strong. 

Expressions  (33)  -  (36)  can  be  manipulated  to  give  the  tempera¬ 
ture,  length,  resistance  and  resistivity  of  the  plasma  armature  in  the 
f o I  low i ng  forms : 
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In  general,  the  R.H.S.  of  these  expressions  are  implicit  functions  of 
the  principal  variables  on  the  L.H.S.  Expressions  (47)  -  (50)  serve  as  a 
self-consistent  set  of  equations  which  must  be  solved  simultaneously  in 
order  to  obtain  the  structure  and  properties  of  the  plasma  armature.  For 
strongly  ionized  plasma  armatures,  the  R.H.S.  is  only  very  weakly 
coupled  to  the  variables  on  the  L.H.S. ,  and  the  above  expressions  show 
in  an  approximate  way  the  explicit  dependence  of  the  plasma  temperature, 
length  and  resistance  on  the  current  and  the  mass  of  the  plasma 
armature,  and  on  the  launcher  parameters  (L’,  M’  and  the  bore  size). 
Indeed,  in  the  case  of  strong  ionization,  we  can  show  from  expressions 
(33)  -  (36)  that  the  following  is  an  approximate  invariant  for  the 
plasma  armature  according  to  our  model: 
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(51) 


where  j  is  the  mean  current  density  in  the  plasma. 


6.  RADIATIVE  HEAT  TRANSFER  TO  THE  WALL 


6.1  INTRODUCTION 

Gross  ablation  by  the  plasma  armatures  of  the  wall  materials 
continues  to  pose  a  serious  problem  for  rail  guns  to  reach  high 
velocity[14]  and  attain  repetitive  firing  capabi I ity [10,13] .  The 
ablation  occurs  even  after  the  arc  has  gathered  speed  from  a  few  hundred 
meters  per  second  to  more  than  a  few  kilometers  per  second. 

The  phenomenon  is  a  result  of  many  complex  processes  occur ing  at 
the  ra i I -armature  and  ra i I - i nsu I ator  interface.  From  classical  arcs 
study,  the  presence  of  cathode  spots  and  sometimes  anode  spots  are  well 
known [41-45] .  In  these  arc  roots,  extremely  high  concentrations  of 
current  occur.  In  many  conventional  applications  of  electric  arcs  in 
which  the  mean  current  density  is  much  smaller  than  that  experienced  in 
railguns,  radiation  intensity  from  the  arc  column  (arc  interior)  is 
relatively  low  and  insignificant.  Under  those  conditions,  the  energetic 
processes  in  the  arc  roots  dominate  the  erosion  of  the  electrodes, 
leading  to  erosion  rates  typically  in  the  range  of  1  -  100  fig  per 
coulomb  of  charge  transfer. 

Unfortunately,  this  does  not  appear  to  be  the  situation  in  the 
case  of  railgun  arcs.  Firstly,  because  of  the  much  higher  current 
density,  the  mean  radiation  intensity  from  the  arc  body  is  relatively 
high.  Secondly,  the  typical  erosion  rates  measured  in  mass  per  unit 
charge  transfer  observed  in  railguns  is  in  the  range  of  10  -  30  mg/C 
[46],  approximately  3  orders  of  magnitude  above  those  observed  in 
conventional  arcs. 


In  this  study  we  pose  a  simple  but  fundamental  question:  Could 
the  mean  thermal  flux  from  the  main  body  of  the  arc  in  railguns  be 
sufficient  to  cause  melting,  vaporization  and  hence  ablation  of  the  wall 
materials?  If  so,  what  are  the  key  parameters  governing  this  mean 
thermal  flux?  Is  there  a  domain  of  suitable  properties  and  structures 
for  the  plasma  armature  so  that  ablation  caused  by  the  mean  thermal  flux 
of  the  arc  could  be  eliminated? 

The  present  work  represents  one  of  a  continuing  series  of 
efforts  to  address  these  issues  in  a  systematic  way.  It  is  a  signifi¬ 
cant  extension  of  an  earlier  study  by  Powell  [7]  in  that  we  explore  the 
answers  to  these  questions  in  a  much  larger  parameter  space.  It  is  also 
extension  of  an  earlier  work  by  Thio[5]  in  that  non-ideal  plasma  effects 
are  included  in  the  present  study. 

The  dominant  contributions  to  the  mean  heat  flux  from  the  arc 
are  radiation,  turbulent  thermal  conduction  and  energy  generated  by 
hyperve I oc i ty  skin  friction  at  the  wall.  It  must  be  noted  that  the 
turbulent  thermal  conduction  and  skin  friction  heat  generation  can  be  as 
large  or  greater  than  the  radiative  transfer  under  certain 
c i rcumstances  [48]  .  Marshal  I  [40]  has  recently  also  drawn  attention  to 
the  presence  of  arc  roots  in  railgun  arcs.  The  result  of  these 
additional  heat  transfer  mechanisms  is  to  aggravate  the  wall  ablation 
problem  in  railguns.  An  excellent  and  parallel  study,  which  carefully 
treats  the  various  forms  of  energy  and  momentum  transfer  to  the  wall 
from  the  arc  through  the  plasma-wall  boundary  layer,  has  been  completed 
recently  by  Tidman,  Goldstein  and  Ainsor[48],  In  this  paper,  we  shall 
focus  our  attention  on  the  radiative  contr i but i on . 

6.2  THE  PHYSICAL  FRAMEWORK  AND  THE  METHOD  OF  ANALYSIS 

Consider  a  stable  plasma  armature  moving  at  constant  accelera¬ 
tion  and  endowed  with  a  uniform  current  density,  a  configuration  of 
minimum  peak  current  concentration  for  a  given  total  current  and  spatial 
extent  of  the  plasma  armature  For  a  gi .en  set  of  gun  parameters  (such 
as  bore  size  and  ra i I  inductance  gradient)  and  plasma  parameters  (such 


as  its  total  mass  and  ionic  species),  using  the  results  developed  in 
Chapter  5,  we  can  calculate  the  temperature  and  the  length  of  the  plasma 
armature.  The  radiative  flux  from  the  arc  on  to  the  rails  and 
insulators  is  then  determined  and  this  in  turn  determines  the  rate  of 
temperature  rise  in  the  materials.  To  this  we  add  the  magnetic  heating 
of  the  rails  caused  by  the  pulsed  current  flowing  in  them.  The  length 
of  the  plasma  armature  determines  its  dwell  time  at  a  given  position  of 
the  wall  as  it  flies  by.  Together,  the  radiative  heat  flux,  the  rate  of 
thermal  energy  generation  by  magnetic  heating,  and  the  dwell  time 
determine  whether  the  melting  point  of  the  material  would  be  exceeded 
during  the  fly-by  of  the  plasma  armature. 

By  appropr i atel y  treating  the  non- idea  I  plasma  behavior  in  our 
model,  the  modelling  results  should  have  a  higher  degree  of 
approximation  to  the  practical  situations. 

6.3  ARMATURE  RESISTANCE,  TEMPERATURE,  AND  LENGTH:  AN  EXAMPLE 

In  this  section,  we  present  the  results  obtained  directly  from 
the  solution  of  the  system  of  equations  (47)  -  (50),  together  with  the 
Saha  equation  (32)  and  the  relevant  supporting  equations  (27)  -  (31), 

(37)  -  (43),  for  an  exemplary  case  of  a  railgun  which  has  a  1-cm  square 
bore  and  an  inductance  gradient  L’  of  0.32  /iH/m.  The  calculations  are 
made  for  currents  between  200  kA  and  400  kA  and  unity  has  been  assumed 
for  all  em i ss i v i t i es .  The  scaling  factor  Sy  has  been  assigned  the  value 
of  0.5  which  is  the  value  used  in  obtaining  the  calibrated  resistivity 
shown  in  Figure  2 . 

Figures  3  and  4  show  respectively  the  estimated  resistance  and 
temperature  of  the  arc  column  versus  the  total  number  of  particles  (i  ns 
and  neutral  atoms)  in  the  armature. 

20 

When  the  number  of  particles  is  small,  say  5  *  10  ,  the 

armature  tends  to  be  short,  the  resistance  tends  to  be  high,  and  the 
temperature  tends  to  be  high. 

In  the  case  of  a  pure  Li  plasma,  the  resistance  is  of  the  order 
of  1  mQ  and  is  relatively  insensitive  to  variation  in  current.  The 
temperature  is  relatively  high-  approximately  25,000  K  at  a  current  of 


L'  -  0.  32  uH/m 
Bore  - 1  cm 


Figure  3 


.06 


.18  .30 


1.2 


Z  4 


Total  Number  of  Particles  NT  (  1022) 


—  Arc  resistance  versus 
ions  and  neutrals)  in 
ra i I  gun  with  an  L ’  of 
in  cross-section. 


the  total  number  of  particles  (charged 
the  plasma  armature  for  the  case  of0a 
0  32  .iH/m  and  a  square  bore  of  1  cm'- 


Total  Number  of  Particles  (  10  / 


Temperature  of  the  plasma  armature  as  a  function  of  tota 
number  of  particles  in  the  armature  for  the  case  of  a 
square-bore  ra!gun  (1  cm  )  #ith  an  L’  of  0  3z  uH/m 


200  kA  and  45,000  K  at  400  kA.  At  these  temperatures,  the  lithium 
plasma  is  essentially  ideal  (7  <<  0.1).  Thus,  the  electron-ion 
contribution  of  the  resistivity  here  could  have  been  well  approximated 
by  Spitzer’s  expression.  Further,  due  to  the  relatively  low  ionization 
potential  of  lithium,  the  Li  plasma  is  nearly  fully  ionized.  Thus  the 
contribution  to  the  resistivity  from  the  electron-neutral  scattering  is 
small  due  to  a  small  concentration  of  the  neutral  particles.  When  these 
conditions  prevail,  the  net  resistance  of  the  armature  behaves  nearly 
independently  of  the  current.  As  the  current  increases,  the  temperature 
increases  and  the  resistivity  drops.  However,  the  length  of  the 
armature  also  decreases  resulting  in  a  smaller  cross-sectional  area 
through  which  the  plasma  current  passes.  The  overall  effect  is  to  make 
the  net  resistance  of  the  armature  insensitive  to  variation  in  the 
current . 

As  the  number  of  particles  in  the  armature  increases,  the  arc 
length  grows,  and  the  resistance  generally  decreases.  In  the  case  of  a 
Li  plasma  armature  in  the  gun  used  in  the  example,  the  armature  resis¬ 
tance  is  reduced  to  a  level  below  0.3  mfi  for  a  current  of  300  kA  when 
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the  number  of  particles  in  the  armature  reaches  a  value  of  2  *  10 
However,  the  plasma  also  cools  with  increasing  number  of  particles.  In 
the  case  of  H  or  0  plasma,  due  to  the  high  ionization  potential  of  the 
atomic  species,  the  degree  of  ionization  in  these  plasmas  decreases 
rapidly  as  the  plasma  cools.  The  consequence  is  the  presence  of  a  high 
concentration  of  neutral  particles  giving  rise  to  a  significant 
contribution  from  the  electron-neutral  scattering  to  the  resistivity  and 
consequently  a  higher  resistance.  Thus  in  these  plasmas,  the  overall 
armature  resistance  at  first  decreases  with  the  number  of  particles 
present  in  the  armature  unti I  a  point  is  reached  when  further  increase 
in  the  number  of  particles  makes  the  neutrals  dominate  the  plasma 
resistivity.  Beyond  that  point,  the  armature  resistance  increases  with 
increasing  number  of  particles. 


For  a  pure  0  plasma  at  300  kA,  the  resistance  decreases  from 

20  22 

approximately  1.2  mQ  for  5  x  10  particles  to  about  0.6  mQ  for  1  x  10 

particles  in  the  armature.  For  a  pure  H  plasma  at  the  same  current,  the 

resistance  appears  to  fal  I  from  1.4  mfl  to  1  mQ  as  the  number  of 

20  22 

particles  is  increased  from  5  x  10  to  1  *  10  When  the  number  of 

22 

particles  increases  beyond  2  /  10  ,  the  armature  resistance  begins  to 

gradually  increase. 

With  reference  to  Figure  4  and  5,  we  see  that  the  armature 
temperature  and  more  importantly  the  radiation  intensity  decreases 
monoton i ca I  I y  with  increasing  number  of  particles. 

20 

At  200  kA  when  the  number  of  particles  increases  from  5  x  10 

22 

to  2  x  10  the  temperature  of  a  Li  plasma  drops  from  25,000  K  to  12,000 

K.  The  corresponding  radiation  intensity  from  the  armature  falls  from 

2  2 

25  MW/cm  to  below  0.1  MW/cm  .  For  an  0  plasma  armature,  the 

temperature  decreases  from  approximately  27,000  K  to  about  13,000  K  and 

2  2 

the  radiation  intensity  goes  from  3  MW/cm  to  about  0.14  MW/cm  .  For  a 
H  plasma  armature,  the  temperature  decreases  from  about  30,000  K  to 

2 

about  20,000  K  and  the  radiation  intensity  drops  from  about  5  MW/cm  to 

2 

approximately  025  MW/cm  . 
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At  400  kA,  when  the  number  of  particles  increases  from  5  *  10 

22 

to  2  *  10  ,  the  temperature  of  a  Li  plasma  drops  from  a  45,000  K  to 

18,000  K.  The  corresponding  radiation  intensity  from  the  armature  falls 

2  2 

from  25  WW/cm  to  below  0.5  MW/cm  .  For  an  0  plasma  armature,  the 

temperautre  decreases  from  approximately  45,000  K  to  about  20,000  K  and 

2  2 

the  radiation  intensity  goes  from  25  MW/cm  to  about  0.8  MW/cm  .  For  a 
H  plasma  armature,  the  temperature  decreases  from  about  47,000  K  to 

2 

about  23,000  K  and  the  radiation  intensity  drops  from  about  30  MW/cm  to 

2 

approximately  1  2  WW/cm  . 

Even  refractory  and  conducting  materials  such  as  tungsten  or 

molybdenum  are  unable  to  survive  the  radiation  from  the  armature  for 

2 

more  than  a  few  us  when  the  radiation  intensity  exceeds  5  WW/cm“  For 

2 

radiation  intensity  below  1  MW/ cm  ,  there  is  a  chance  for  these 
materials  to  survive  if  the  dwell  time  of  the  plasma  armature  is  of  the 


Radiation  Intensity  (MW/cm  ) 


order  of  tens  of  fis .  The  thermal  conductivity  of  the  materials  plays  an 
important  part.  In  the  next  section,  the  radiative  heating  of  the  wall 
materials  is  analysed. 

6.4  PLASMA  RADIATIVE  HEATING  OF  THE  WALL  MATERIALS 

The  wall  of  the  bore  is  subjected  to  an  incident  radiative  flux 
g i ven  by , 

Q  =  e  (52) 

where  T^  is  the  arc  temperature,  <7^  is  the  Stefan  constant,  and  e  the 
effective  emissivity.  We  shall  assume  unity  for  the  effective 
emissivity  for  want  of  exact  experimental  data.  This  will  tend  to  over¬ 
estimate  the  radiative  flux.  The  temperature  T^  strictly  should  be  the 
temperature  of  the  arc  at  a  point  which  is  approximately  one  optical 
(Rosseland)  path  length  from  the  wall.  Again,  we  can  obtain  a  worst- 
case  estimate  by  using  the  arc  column  temperature  calculated  from  the 

model  as  the  value  for  T. . 

A 

Before  the  onset  of  melting,  the  wall  temperature  increases 
according  to  the  standard  equation  for  heat  conduction  in  solid, 


where  0  is  the  mass  density,  c  is  the  specific  heat  per  unit  mass  of 
the  rail,  *  is  its  thermal  conductivity,  and  x  is  the  displacement  into 
the  wall  perpendicular  to  its  surface.  The  boundary  condition  at  the 
wall  surface  (x  =  0)  providing  the  radiative  heat  flux  is, 

T*  =  -  Q  =  -  e  csj4  <54) 

Equation  (53)  with  the  boundary  condition  (54)  can  be  solved  to  give  the 
temperature  at  a  point  x  and  time  t  to  be, 


t>-S 


(55) 


r(x,t)  =  T  ♦  Q  {(ttX)  exp  (“X  x  )  -  x  Erfc  (Ax)}/* 

x  =  (/?  Cy/^t)1^2 


where  Tg  is  the  initial  temperature  of  the  ra i 
x  =  0,  the  temperature  at  time  t  is. 


At  the  ra i I  surface , 


T(0,t)  =  TQ+  2  Q  (t/xpc^K)1/2 ,  t<tr 


In  the  above,  the  time  t  is  measured  from  the  instant  the  plasma 
armature  first  reaches  the  point  of  interest.  The  above  expression 
applies  only  during  the  fly-by  of  the  plasma  armature.  The  dwell  time 
of  the  plasma  is  indicated  as  tr .  After  this  time,  the  wall  temperature 
would  then  fall  off.  It  can  be  shown  that,  after  the  armature  fly-by, 
the  surface  temperature  of  the  wall  at  time  t  >  tr  is  given  by 


T  (0 ,  t)  =  Tg+  2  Q  {t1/2-  (t  -  tr)1/2}/(^cv*)1/2 


The  above  heating  of  the  wall  by  the  plasma  radiation  applies  to  both 
the  rail  and  the  insulator.  For  the  rail,  an  additional  heating 
mechanism  needs  to  be  included:  the  magnetic  heating  (or  commonly 
called  Joule  heating)  of  the  rails.  The  magnetic  and  radiative  heating 
of  the  rails  strictly  are  a  coupled  thermal  process.  In  a  recent  paper, 
Powell [7]  treats  this  coupled  phenomenon  by  carefully  formulating  a  1-D 
model  in  which  the  equations  governing  the  velocity-skin  effect 
(magnetic)  heating  of  the  rail  is  solved  simultaneously  with  the 
equations  for  the  radiative  transfer.  As  also  noted  by  Powell,  there 
exists  at  present  a  high  degree  of  uncertainty  between  application  of 
the  theoretical  results  and  interpretations  of  experimental  observations 
in  any  exact  form.  In  view  of  this  uncertainty,  it  appears  that  a 
consistent  degree  of  accuracy  can  still  be  obtained  even  if  we  seek  a 
simplification  of  the  mathematical  problem  by  un-coupling  the  magnetic 


m 


heating  of  the  rail  from  the  radiative  heating.  In  this  procedure,  we 
adopt  a  worst-case  approach  in  estimating  the  heating  effects  from  each 
of  these  contributions.  By  treating  the  combined  effects  as  a 
superposition  of  the  separate  effects,  the  overall  estimate  for  the 
temperature  rise  in  the  rails  can  reasonably  be  expected  to  be 
conservative  in  the  maximal  sense.  This  simplification  allows  rapid 
parametric  studies  over  a  wide  range  of  plasma  and  rail  parameters  to  be 
performed.  The  approach  is  in  essence  equivalent  to  a  I ump-parameter 
method . 

6.5  MAGNETIC  (OR  JOULE)  HEATING  OF  THE  RAILS 

The  precise  way  of  taking  account  of  magnetic  heating  of  the 
rails  is  to  solve  the  coupled  set  of  partial  differential  equations 
governing  strong  magneto-thermal  diffusion  in  a  conductor , coup  I ed  to  the 
electrodynamics  of  the  moving  arc.  An  upper  bound  to  the  temperature 
rise  due  to  magnetic  heating,  however,  can  be  obtained  as  follows.  The 
magnetic  field  within  the  plasma  armature  decreases  from  the  breech-end 
of  the  armature  to  the  muzzle-end  of  the  armature  according  to  the 
uniform  current  density  model  developed  above.  During  the  fly-by  of  the 
plasma  armature,  at  a  given  station  on  the  surface  of  the  rail,  one 
would  see  a  magnetic  field  increasing  monoton i ca I  I y  with  time.  In  the 
case  of  an  uniform  current  density  model,  this  monotonicity  degenerates 
into  one  of  linearity.  After  the  plasma  armature  has  completely  passed 
by,  the  magnetic  field  takes  on  the  value  of  the  field  produced  by  the 
constant  current  flowing  in  the  rails.  The  temperature  rise  due  to  this 
temporal  monotonic  increasing  variation  of  magnetic  field  with  a  plateau 
amplitude  after  a  certain  lapse  of  time  is  always  less  than  that  due  to 
a  sudden  switch-on  of  a  magnetic  field  with  the  same  plateau  amplitude. 
For  a  conservative  (maximal)  estimate  of  the  rail  surface  temperature 
due  to  Joule  heating,  we  may  therefore  use  the  result  for  the  case  of  a 
sudden  switch-on  of  a  magnetic  field. 

For  the  sudden  switch-on  of  a  magnetic  field,  the  peak 
temperature  rise  occurs  at  the  conductor  surface  and  at  the  instant  of 


switching  on  the  field,  if  effects  associated  with  temperature-dependent 
material  properties  were  neglected. 

The  question  remains  whether  2-D  or  even  3-D  effects  are 
important  in  determining  the  magneto-diffusion  heating  effects  in  the 
case  of  finite  and  curved. ra i I s .  It  can  be  shown  that,  in  a 
sufficiently  short  time  scale  and  in  a  sufficiently  thin  layer  of  the 
rail  closed  to  the  surface,  2-D  or  3-D  magneto-diffusion  may  be  ignored. 
In  particular,  the  peak  surface  temperature  which  occurs  at  t  =  0  may  be 
obtained  from  a  solution  of  a  1-D  magneto-diffusion  problem: 


with  the  appropriate  boundary  conditions,  where  /i  is  the  magnetic  perme¬ 
ability,  rj  is  the  rail  resistivity  and  the  other  symbols  have  their 
usual  meaning.  Again,  locally  and  for  the  very  short  time  of  interest, 
the  boundary  conditions  can  be  chosen  to  be  the  same  as  those  for  a 
semi  - i nf i n i te  conducting  half-space,  with  the  surface  of  the  half-space 
corresponding  to  the  rail  surface.  The  problem  now  reduces  to  one 
originally  solved  by  Kidder[49]  and  analyzed  by  several  other  investi¬ 
gators,  e.g.[47].  From  their  solution,  the  maximum  surface  temperature 
rise  can  be  given  as, 
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(59) 


We  use  this  expression  to  provide  the  upper  bound  to  the  magnetic 
heating  of  the  rail  due  to  the  moving  plasma  armature. 


6.6  MINIMUM  ARMATURE  VELOCITY  FOR  WALL  SURVIVABILITY 


If  the  plasma  armature  enters  a  barrel  section  with  an  entry 
velocity  Vg  and  is  subjected  to  a  constant  acceleration  a,  the  dwell 
time  t  of  the  armature  at  the  entry  point  is  given  by, 

tr  =  {(vQ2  ♦  2a  I ) 1/2  -  v0}/a  (60) 

where  I  is  the  length  of  the  plasma  armature  which  is  obtained  from 
solving  the  plasma  model  as  described  by  Chapter  5.  For  a  constant 
acceleration,  the  armature  dwell  time  decreases  with  increasing  entry 
velocity  Vq.  In  the  case  of  the  rail,  the  total  temperature  rise  is 
given  by, 


AT 


ra  i  I 


(61) 


where  T  is  the  temperature  rise  due  to  radiative  heating  given  by 
expression  (56),  ant  T(  is  the  temperature  rise  due  to  Joule 
(magnetic)  heating  given  by  expression  (59).  In  the  case  of  insulator, 
the  total  temperature  rise  consists  of  only  the  radiative  term.  From 
expression  (56),  we  see  that  the  total  temperature  rise  of  the  wall 
decreases  with  decreasing  armature  dwell  time.  In  order  to  control  the 
wall  temperature  rise  to  be  below  a  certain  specified  value,  the 
armature  dwell  time  is  required  to  be  below  a  certain  value.  In  turn, 
this  requires  a  certain  entry  velocity  for  a  given  acceleration  of  the 
armature.  Substituting  expression  (60)  for  the  dwell  time  in  the 
radiative  heating  expression  (56),  and  introducing  the  resultant 
expression  in  (61)  and  using  the  Kidder’s  temperature  expression  (59), 
if  necessary,  in  expression  (61),  an  equation  is  obtained  for  the 
minimum  entry  velocity  for  a  given  acceleration,  or  a  minimum 
acceleration  for  a  given  entry  velocity. 

Let  T  be  the  meftinq  point  of  the  wa I  I  material  and  Tn  its 
m  3  0 

initial  temperature.  The  maximum  allowable  temperature  rise  at  the  wall 
surface  before  onset  of  melting  is  given  by  AT  =  T  -  T-. .  In  practice, 


because  of  thermal  softening  of  the  wall  material  at  elevated  tempera¬ 
ture  and  the  requirement  to  withstand  the  high  pressure  of  the  plasma 
armature,  the  allowable  temperature  rise  for  the  wall  is  considerably 
lower  than  is  given  by  the  above  expression  (based  upon  the  melting 
point)  in  order  to  avoid  significant  damage  of  the  wall  material. 

Figure  6  shows  the  result  of  solving  the  resultant  equation  for 
minimum  entry  velocity  to  control  the  temperature  rise  in  the  tungsten 
rails  to  below  2500  K  for  a  railgun  with  1-cm  square  bore,  0.32  ^H/m, 
and  a  1-g  projectile.  The  calculations  are  made  for  the  3  different 
single-element  plasmas,  Li,  0  and  H,  and  for  different  values  of 
current,  200  kA  and  300  kA.  The  minimum  entry  velocity  is  shown  versus 
the  total  number  of  particles  in  the  armature. 

At  300  kA  for  this  gun,  the  required  velocity  for  the  armature 

to  avoid  rail  damage  can  be  seen  to  be  very  high  (exceeding  50  km/s)  for 

21 

all  the  3  different  plasmas  when  the  armature  contains  less  than  10 

22 

particles.  Only  when  the  armature  contains  more  than  10  particles 

does  the  minimum  armature  velocity  fall  below  20  km/s  in  the  case  of  a  H 

plasma  armature,  and  below  10  km/s  in  the  case  of  an  0  plasma.  In  the 

case  of  a  Li  plasma,  the  calculation  indicates  that  an  entry  velocity  as 

low  as  1  km/s  would  be  sufficient  to  avoid  radiative  damage  if  more  than 
22 

10  particles  are  present  in  the  plasma  armature. 

22 

At  200  kA,  the  situation  improves  markedly.  With  more  than  10 

particles  in  the  armature,  the  tungsten  rails  appear  to  be  able  to 

survive  the  armature  radiation  for  moderate  armature  velocity  for  all 

the  3  different  plasmas.  Specifically,  for  a  H  plasma  a  minimum 

armature  velocity  of  approximately  5  km/s  is  required,  and  for  both  0 

and  Li  plasmas,  an  entry  velocity  of  the  order  of  100  m/s  appears  to  be 

all  that  is  neccessary  to  avoid  radiative  damage.  On  the  other  hand,  if 
21 

only  10  particles  are  used,  entry  velocity  of  8  km/s,  10  km/s,  and 
20  km/s  are  required  for  a  Li,  0  and  H  plasma  armature  respectively. 

6.7  CONCLUSION 

On  the  one  hand,  the  modelling  results  show  that  arc  radiation 


is  sufficiently  high  to  cause  ablation  under  Lne  conditions  of  most 


railgun  experiments  conducted  in  the  past.  Under  the  same  conditions 
the  arc  also  tends  to  have  a  rather  high  resistance  (above  1  mfi) .  These 
properties  were  a  consequence  of  operating  the  plasma  armature  under 
less  than  optimum  conditions.  In  particular,  the  plasma  armatures  in 
these  experiments  were  initiated  with  less  than  the  required  number  of 
particles  to  prevent  radiative  ablation.  As  a  result,  the  initial 
plasma  is  very  hot,  emits  radiation  at  a  high  intensity,  and  begins  to 
ablate  the  wall  materials.  Also,  because  of  insufficient  number  of 
particles,  the  arc  has  a  short  length  giving  rise  to  a  high  resistance. 
On  the  ether  hand,  the  modelling  results  indicate  that  there  indeed 
exist  domains  of  operating  parameters  in  which  radiative  ablation  of  the 
wall  materials  can  possibly  be  eliminated.  For  a  given  set  of  gun 
parmeters  and  a  given  current,  we  can  generally  find  a  combination  of 
number  of  particles  in  the  armature  and  the  armature  velocity  by  which 
radiative  ablation  would  not  occur.  We  find,  indeed,  that  by  using  a 
low  atomic  weight,  low  ionization  potential  ionic  species  such  as  Li, 
the  armature  velocity  could  be  as  low  as  a  few  km/s  when  radiative 
ablation  can  be  totally  avoided  under  typical  railgun  operating 
conditions,  and  that  the  mass  of  the  lithium  armature  required  to 
initiate  the  arc  remains  below  10??  of  the  projectile  mass.  In  general, 
a  large  number  of  particles  are  required  to  produce  a  sufficiently  cool 
plasma.  The  greater  length  of  the  resulting  armature  in  general  also 
tends  to  lower  its  electrical  resistance. 

Whether  such  armatures  can  be  conveniently  produced  and 
maintained  throughout  the  launch  remains  a  development  for  the  future. 
Further,  the  greater  length  for  the  armature  could  possibly  increase  the 
overall  viscous  drag,  and  the  lower  conductivity  of  the  cool  plasma 
could  possibly  lead  to  weaker  coupling  between  the  magnetic  field  and 
the  plasma.  Both  effects  could  contribute  to  its  instability. 
Preliminary  investigations  on  these  issues  have  been  made  on  the  related 
program[16] .  Much  more  theoretical  and  experimental  research  on  these 
issues  are  needed. 
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